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FREQUENCY-DOMAIN  UPPER  BOUNDS 
FOR  SIGNALS  ON  A  MULTICONDUCTOR  TRANSMISSION  LINE 
BEHIND  AN  APERTURE 

WILLIAM  A.  DAVIS  and  MOHAMMAD  K.  SISTANIZADEH 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 

Abstract 

This  paper  develops  a  technique  for  bounding  the  maximum  voltages 
and  currents  at  terminations  of  a  multiconductor  transmission  line  (MTL) 
located  behind  an  aperture-perforated  conducting  screen  excited  by  an 
electromagnetic  field  in  the  frequency  domain.  The  electromagnetic  field 
is  coupled  through  a  small  aperture  as  the  excitation  of  a  multiconductor 
transmission  line  behind  the  aperture.  A  model  is  presented  in  terms  of 
external  and  internal  sources  which  in  turn  create  traveling  waves  on  the 
multiconductor  transmission  line.  These  traveling  waves  transfer  energy 


to  the  terminations.  The  energy  at  a  termination  is  translated  into 
voltages  and  currents  from  which  the  upper  bounds  are  determined.  These 
upper  bounds  are  obtained  using  vector  norms  and  associated  matrix  norms. 
The  formulation  is  presented  in  the  frequency  domain  to  obtain  useful 


upper  bounds  for  analysis  of  multiconductor  transmission  line  geometries 


with  aperture  excitation. 


Introduction 


In  designing  some  systems,  the  designer  should  be  able  to  character¬ 
ize  the  penetration  of  electromagnetic  pulses  (EMP)  or  lightning  signals 
through  apertures  of  general  shapes  as  well  as  quantify  the  effects  of 
the  coupled  energy  on  transmission  lines  located  in  the  vicinity  of  the 
aperture.  Apertures  that  are  of  concern  to  the  designer  are  usually 
electromagnetically  small  over  the  spectrum  of  the  IMP,  or  lightning, 
and  their  existence  may  be  for  some  purpose,  e.g.,  windows,  open  access 
holes,  or  they  may  be  unintentional  as  in  the  case  of  cracks  around 
doors  or  plates  covering  access  ports  or  poor  electrical  seams.  The 
analysis  of  the  coupling  (penetration)  problem  has  been  investigated 
by  a  large  number  of  people  since  1897.  The  first  scientist  to  propose 

a  solution  was  Lord  Rayleigh  [1],  whose  solution  was  expressed  as  an 

2n 

ascending  power  series  of  the  wavenumber  k(  *  —  )  where  X  is  the  wave¬ 
length.  Others  include  Bethe  [2],  Bouwkamp  [3],  and  more  recently 
Butler  [4]. 

The  coupling  of  the  energy  from  an  incident  electromagnetic  wave 
to  a  transmission  line  located  behind  an  aperture-perforated  conducting 
screen  has  also  been  investigated  by  many  engineers  and  physicists  in 
the  past  decade.  Kajfez  [5]  has  computed  the  coupled  energy  by  the  use 
of  equivalent  electric  and  magnetic  dipole  moments  along  with  both  mode- 
matching  and  reciprocity  techniques  to  obtain  equivalent  sources  of  a 
transmission  line  model.  Butler  and  Umashankar  [6]  have  approached  the 
problem  numerically  by  the  method  of  moments,  and  have  formulated  lntegro- 
differentlal  equations  for  a  finite-length  wire  with  arbitrary  orientation 
behind  an  arbitrarily  shaped  aperture.  Davis  [7]  has  developed  a  method 
for  bounding  the  maximum  voltage  and  current  levels  at  terminations  of  a 


single  wire  behind  an  aperture. 

This  paper  extends  the  bounding  problem  of  a  single  wire  to  the 
problem  of  obtaining  an  upper  bound  for  the  computation  of  the  voltages 
and  currents  at  terminations  of  multiconductor  transmission  lines  (MTL) 
located  behind  an  aperture-perforated  conducting  screen.  A  brief  summary 
of  the  system  model  will  be  given,  followed  by  a  more  detailed  discussion 
of  the  signal  computations  and  the  resultant  upper  bounds. 


System  Model 

A  typical  multiconductor  transmission  line  (MTL)  problem  of  interest 
is  shown  in  Fig.  1.  The  aperture  A  is  the  source  of  coupled  energy 
which  is  modeled  by  a  pair  of  electric  and  magnetic  current  dipoles  above 
a  closed  aperture  as  shown  in  Fig.  2.  It  is  easily  shown,  from  the  small 
aperture  theory  as  used  by  Kajfez  [8],  that  these  point  source  amplitudes 
are  given  by 


—  A  g  Q—  A 

c  *  c  y  ■  iu>ea  E  y 

e  eyJ  J  e  y  7 


a) 


and 


a  SC  — A 

c  *  c  x  -  -iojya  H  x 
m  mx  J  m,xx  x 


(2) 


respectively,  where  at  is  radian  frequency,  a  and  a  are  the  required 

6  m }  xx 

polarizabilities  of  the  aperture  and  Esc~  and  Hsc~  are  the  electromagnetic 

sc— 

fields  below  the  aperture  with  no  aperture  present.  The  is  of  no 

importance  to  the  problem  and  has  been  neglected^ as  have  the  low  level 

sc—  sc— 

aperture  fields  of  the  MTL.  In  the  process  of  bounding,  E  and  H 

y  x 

are  assumed  to  be  given  while  the  a&  and  ^  may  either  be  given  by  the 
appropriate  geometry  related  coefficients  or  bounded  by  the  polarizabi¬ 
lities  of  an  ellipse  which  circumscribes  the  aperture  [7]. 

Ue  will  restrict  the  MTL  to  be  in  a  vacuum  environment  and  to 


consist  of  N  parallel  wires.  In  such  a  case  there  are  N  transverse 


Figure  1 


A  mulciwire  transmission  line  (MTL)  parallel 
to  a  plane  with  aperture  A. 


Figure  2.  Model  of  the  aperture  as  replaced  by 
equivalent  current  dipoles. 


electromagnetic  (TIM)  modes  possible  on  chs  MTL.  The  equations  that 

govern  these  nodes  are 


£  |v  >  -  -JwL| I  > 


(3) 


and 

£  1 1  >  -  ->K|V  >  (4) 

vhtrt  |v  >  and  |l  >  ara  vactors  raprasanting  the  N  wire  voltages  and 
N  z-directed  wire  currants.  The  induction  matrix  L  and  capacitive  co¬ 
efficient  matrix  K  represent  respectively  the  inductive  and  the  capacitive 
affects  of  the  MTL,  The  current  | I  >  is  equal  to  (cjQL>  )  for  z-directed 
propagation,  c  *  speed  of  light,  which  gives  the  corresponding  charge 
definition 

|ql  >  -  k|v  >  (5) 


Thus  K  may  be  found  from  a  two-dimensional  capacitive  boundary  value 
problem.  The  matrix  L  is  then  given  by  the  Inverse 


L-^R*1  .  (6) 

c 

The  solutions  to  (3)  and  (4)  may  be  written  as 
N 

|V  >  -  l  (a.e'j8z  +  b.eJ62H.  >  (7) 
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and 
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i-1 

where  S  «  <*>/c  and  the  summation  is  over  the  N  voltage  modes  \<p^  >.  The 
>  ara  orthogonal  and  are  normalized  in  the  following  inner  product 


sense  for  unit  power: 


>  -  «ij. 


(9) 


It  is  convenient  to  choose  the  |<t>^  >  to  be  the  eigenvectors  of  K  to  give 
XjC  <  ♦jUj  >  -  5±j  (10) 

where  is  the  eigenvalue  of  K.  With  these  definitions,  the  net 
power  at  any  point  on  the  MTL  is  given  by 

p  -  *S  !  (|aj2  -  |b  I2).  (11) 

i-1  1  1 

Eqs.  (7)  and  (8)  may  be  written  in  more  compact  form  as 

|v  >  -  My(|a  >  e"j6z  +  |b  >  eJ8z)  (12) 

and 

I1  >  -  cKHy(|a  >  e‘j6z  -  |b  >  ej8z)  (13) 

and  modeled  as  in  Fig,  3  where 

My  *  [I*!  > . UN  >  ]• 

The  matrices  and  of  Fig,  3,  represent  the  phase  delay  between  the 

aperture  and  terminations,  JT^  and  r^.  The  traveling  wave  sources  |ag> 

and  |b  >  represent  the  aperture  coupling  to  the  MTL. 
s 

The  sources  are  obtained  using  the  reciprocity  relation 
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where  the  m  subscript  refers  to  magnetic  current  sources.  For  each 
voltage  mode  >  the  corresponding  electromagnetic  field  intensities 
are  e^  and  h^  respectively,  where  for  n  ■  /y?e 


■  V  f  i  Ty~^  ii 


The  total  I  and  H  have  expansions  similar  to  (7)  and  (8).  Letting  E  ,  H 

d  i 

J  ,  and  J  be  equal  to  E,  H,  c  6(r),  c  6(r)  respectively  (6(r)  the  3-D 
A  ma  6  m 

Dirac  delta),  then  using  e^  and  with  e"^z  for 

aSi  -  ‘st-Cey  -  (16) 


h  and  “fa* 


and 
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or  equivalently 
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(17) 


(18) 


and 
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(19) 


Once  the  distances  and  and  the  termination  voltage  reflection 
matrices  T_^  and  are  known,  the  total  termination  voltages  may  be 
obtained.  These  reflection  matrices  are  given  by 


and 


|b  >  e+jg*4  -  rja>e'j3*4 


|a>e^^3,r3|b>e+J^3 


and  the  delay  matrices  by 


4  -  e'^4  I 


and 


(20a) 


(20b) 
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where  I  is  the  N  *  N  identity  matrix.  Concentrating  on  the  termination 


x  I|ag  >  +  I3e~:i26i!'3|bs>  ]e"jSJl4.  (22) 

is  the  matrix  of  (12)  and  the  second  term  accounts  for  the  sum  of 
| a  >  and  |b  >  .  The  inverse  term  accounts  for  the  multiple  reflections 
cn  the  line  with  l  corresponding  to  the  total  line  length.  The  remaining 
terms  account  for  the  source.  The  corresponding  current  may  be  obtained 
simply  by  replacing  +  T^]  of  (22)  by  cKMy[I  -  T^]. 


Bounds 

The  problem  may  now  be  stated  as  obtaining  the  upper  bound  of  the 
voltage  AV  between  any  two  wires  at  z  *  JL^.  This  voltage  can  not  exceed 
twice  the  voltage  of  any  wire  to  the  ground  plane.  Thus 


AV 


<  2 
max  — 


'V 


where  the  infinity  norm  of  |V^  >  is  the  maximum  absolute  value  (V^l  over 
the  wires.  However,  the  infinity  norm  of  a  vector  is  less  than  or  equal 
to  the  two  norm  given  by 


IV 


/|V4ll 


|V4„I 


to  give 


AV  <  2 
1  ‘max  — 


(23) 


To  complete  the  bounding  process,  three  important  properties  of  matrix 


norms  which  follow  are  needed: 


I!  A  ||  *  0  if  and  only  if  A  =  0 
II  A  +  B  ||  <  ||  A  ||  +  ||  B  || 

II  A  B  ||  <  || A||  ||  B  ||  , 

along  with  the  two  norm  of  a  matrix 

||  A  ||  2  31  [max (eigenvalue  A^A]^ 

where  A"^  is  the  conjugate  transpose  of  A. 

For  passive  terminations,  the  norms  of  both  and  _T^  are  less  than 
or  equal  to  unity.  With  the  norm  of  given  in  terms  of  A^  as 

||M„  ||  =  l//c  min(A  .)  , 

2  i 

the  upper  bound  on  AV  becomes 


AV  <  - 
'max 


a  > 
1  s 


b  > 


m^nC^) 


(1  -  e"°T) 


(24) 


where  represents  the  minimum  round-trip  loss  on  the  MTL  associated  with 
the  multiple  reflections.  Substituting  for  |as>  and  |bg  >  from  (18),  (19), 
(1),  and  (2),  this  upper  bound  becomes 


AV  < 

1  'max  — 
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The  corresponding  maximum  current  1 1 1 mav  is  simply  given  by 


c  AV 


'max 


/max(A  )min(A  )  . 
i  i 


max 
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(26) 


As  a  simple  example  of  using  (25) ,  we  will  consider  the  dual  strip¬ 
line  geometry  of  Fig.  4  as  described  by  Kajfez.  The  K  for  this  geometry 


may  easily  be  obtained  from  two-dimensional  electrostatics  [8]  as 


£  =  eo 


15 

-10 

-10 

15 

with 


^1  3  5e0  *  *2  *  25eo  * 


The  appropriately  normalized  voltage  eigenvectors  of  K  are  given  by 


*1*  -  /lo 


/  50 


1 

-1 


The  normalized  electric  field  term  is  obtained  by  multiplying  the  |<j>i> 

by  K  to  obtain  the  corresponding  (Q  >  from  (5).  Reverting  to  the 

Li 

electrostatic  analysis,  the  short  circuited  field  at  the  aperture  may  be 
obtained  to  give 


■V  3  VoT 


12 

>22 


-307.0 

137.3 


The  remaining  parameters  are  chosen  to  be 

.-3 


O'p 

r»SC“ 


nHSC* 


90.652  ^  10  (or  ldB/100  ft.) 
100  kV/m 


.sc- 


and 

a  ■  aperture  radius  ■  1  cm 

where  *  120ir  and  the  line  lengths  and  X ^  are  7  and  5  meters 
respectively  as  used  by  Kajfez  [8].  A  circular  aperture  of  radius  1  cm 


results  in  polarizabilities  of  o  -6.67  *  10  ^  and  a 


m,xx 


*  2a  [7].  The 


short  circuited  fields,  ESC  and  Hsc  ,  correspond  to  a  plane  wave 

y  ^ 

traveling  along  the  surface. 

Substituting  in  (25),  the  termination  voltage  for  a  radian  frequency 
cj  will  not  exceed 


|AV|  £  (2.38  x  10'7)o)  V. 


This  forms  a  useful  frequency-domain  bound  for  the  problem  presented  with 
a  voltage  less  than  10  Volts  for  frequencies  below  6.7  MHz. 

For  comparison,  this  problem  has  been  solved  exactly  for  open- 
circuit  terminations  on  the  MTL.  For  such  a  case 


Determining  the  |a  >  and  |b  >  of  Eq.  (22)  from  (1),  (2),  (18)  and  (19), 
s  s 

the  results  were  computed  and  are  plotted  in  Fig.  4  along  with  the  bound. 
A  modified  bound  is  also  plotted  which  represents  the  actual  bound  of 
the  particular  problem.  The  difference  in  bounds  is  4.1  which  seems 
slightly  unreasonable  until  the  bounding  approach  is  examined.  A  factor 
of  two  arises  in  the  bound  to  account  for  a  differential  mode  which  does 
not  occur  in  the  case  considered.  The  triangle  inequality  used  in  the 
bound  of  | a  >  and  |b  >  contributes  another  1.5.  The  product  of  the 
2-norms  of  and  |e^>  versus  the  »-norm  of  (Myle^)  contributes  a  1.29 

factor.  A  small  contribution  also  occurs  due  to  some  of  the  neglected 
loss  terms.  In  light  of  these  observations,  the  resultant  bound  is  very 


reasonable 


O'*  w  Sec./ ^4^) 


Conclusion 


This  paper  has  presented  an  approach  to  bounding  signal  levels  at 
MTL  terminations  due  to  excitation  by  an  aperture.  The  approach  requires 
knowledge  of  the  external  fields,  maximum  aperture  dimension,  cross 
section  of  the  MTL,  and  the  length  of  the  MTL.  The  results  use  no  know¬ 
ledge  of  the  terminations.  Such  an  approach  is  useful  for  system  hardness 
evaluations  in  aircraft  and  other  systems  to  BMP  or  other  incident  energy. 
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